Orbital effect of in-plane magnetic field on quantum transport in chaotic lateral dots 
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We show how the in-plane magnetic field, which breaks 
time-reversal and rotational symmetries of the orbital mo- 
tion of electrons in a heterostructure due to the momentum- 
dependent inter-subband mixing, affects weak localisation 
correction to conductance of a large-area chaotic lateral quan- 
tum dot and parameteric dependences of universal conduc- 
tance fluctuations in it. 



A high sensitivity of phase-coherent transport through 
quantum dots to external perturbations has recently en- 
abled one to transform studies of mesoscopic effects |^-^] 
into a spectroscopic tool for detecting tiny energetic 
changes in the electron gas and for studying electron 
dephasing and inelastic relaxation rates A conve- 

nient object 01, once used as a mesoscopic thermome- 
ter [H, consists of a lateral semiconductor dot weakly 
coupled to the reservoirs via two leads, I and r, each 
with Ni_r ^ 1 open conducting channels, and, there- 

fore, quantum conductances gi,r — —Ni^r- Informa- 
tion concerning fine energetic characteristics of single- 
particle electron states in a dot can be extracted from 
the variance and parametric correlations of universal con- 
ductance fluctuations (UCF), Sg = g — {(g)), measured 
as random oscillations of the dot conductance, g around 
the mean value, ((g)) = gWr/igi + upon variation 
of a perpendicular magnetic field [^jj^ , the Fermi energy 
|,| or the dot shape 0. 

Energetic resolution of such a spectroscopy is set by the 
level broadening of single-particle states in a particular 
device, which is limited by the carrier escape into the 
leads, 

= [Ni + Nr)A/h, 

where A = 2Trh? /mS is the mean level spacing of single- 
particle states of spin-polarized electrons with mass m in 
a dot with area S. The use of larger dots with weaker 
coupling to the leads increases the sensitivity of the dot 
conductance to the variation of external parameters. The 
use of larger dots also enables one to assess directly the 
low excitation energy characteristics of the 2D electron 
gas, since the electron properties in l-^lO/xm^-area dots 
containing 10'^ — 10* particles are less affected by the con- 
finement effects. Recently, large area dots were used for 
studying spin-polarisation of a 2D electron gas [Q. In 
order to enhance coupling between a magnetic field and 



electron spin, J. Folk et al Q used a magnetic field finely 
tuned to be orientated exactly parallel to the plane of 
2D electrons. One observed a strong suppression of the 
variance ((i5f;^(-B|| ))) by an in-plane field interpreted 
in terms of a spin-orbit coupling in the 2D gas ||^ . This 
feature has also been accompanied by an observation px| ] 
of such a negative weak localisation (WL) magnetore- 
sistance caused by an in-plane field that one would 
relate to the time-reversal symmetry breaking in the or- 
bital motion of electrons. In the present publication, we 
assess to which extent one can reduce the infiuence of an 
ideally in-plane tuned magnetic field on quantum trans- 
port in lateral semiconductor dots to spin effects alone, 
that is, we determine the range of fields -B|| that would 
affect WL and UCF's in experimentally studied devices 
Q , in addition to spin-related phenomena. 

The influence of an in-plane magnetic fleld on the or- 
bital motion of carriers in a heterostructure or quantum 
well is the result of a flnite width, Xz of a 2D layer, and it 
has been discussed previously in various contexts ||ll|,|l2| . 
The Lorentz force generated by a planar fleld on electrons 
moving across within the 2D plane mixes up the elec- 
tron motion along and across the conflncment direction, 
thus resulting in the electron momentum, p dependent 
subband mixing and, therefore, in a modiflcation of the 
2D dispersion, E{p) — ^ E{B\^,p). In particular, (a) the 
2D electron mass increases in the direction perpendicular 
to , whereas (b) in heterostructures which have no in- 
version symmetry in the form of confining potential, _B|| 
also lifts the p symmetry in the dispersion law [ pd] |: 

i?(S||,p) - E{B\i,~p) cx (p. [4 X ^ 0. 

The latter change in dispersion has potential to reduce 
the fundamental symmetry of chaotic dot from orthogo- 
nal (o) to unitary (u), as a perpendicular magnetic field 
would do. Below, we show that the efflciency of time- 
reversal symmetry breaking by an in-plane magnetic field 
can be characterized using the rate t^|| 5i3y -I- aB| de- 
scribed in Eq. (^. A similar conclusion has recently 
been made in Ref. jl^]. Without spin-orbit effects, this 
parameter would determine the value of the WL correc- 
tion, gwL(^ii) = {{9{B\\))) — {{g))u and of the variance of 
UCF, (((55^(i3||))), as compared to their nominal values. 



gwL(O) ^ {{g))o - {{9))u and {{Sg^))u MM'- 



5wl(S||) = 5wl(0) 1 + r^ij/r^sc 



{{Sg'{Bn))) = {{Sg'))Jl+ l + r„!/r_i 
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The latter parameters can be studied from the UCF's fin- 
gerprints measured in the 'shape of a dot' space in multi- 
gate devices , or by varying the Fermi energy in back- 
gated dots. The rise in the mass anisotropy upon the in- 
crease of would also manifest itself: as a change in the 
UCF pattern. A varying dispersion relation for electrons 
studied at different fields, and B||2i can be charac- 
terized using the rate r^^(i3||i, B||2) oc (Bp — -S|2)^ 
Eq.(^, which can be used to describe auto-correlation 
properties of a full B^-dependent UCF pattern, 



(2) 



In the presence of an in-plane magnetic field, the ef- 
fective 2D Hamiltonian for electrons in a heterostructure 
with a potential profile V{z) can be obtained from the 
3D Hamiltonian, 

H,„.-gi.n=)+ '-"-^"'% »K»). (3) 
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using the plane wave representation, "iip — e^^^^^ipQp{z) 
for the lowest subband electrons. Here, A = {z — Zi))B\\ x 
Iz is the vector potential, zq = (0|z|0) is the center of 
mass position of the electron wave function |0) = tf^^\z) 
in the lowest subband for B\\ — 0, and u{r, z) is a combi- 
nation of Coulomb potential of impurities and lateral po- 
tential forming the quantum dot. Due to mixing between 
subbands |0) and \n > 0) by an in-plane magnetic field, 
z-dependent components (pop(z) are different for different 
in-plane momenta, p, and we use both the perturbation 
theory analysis and a numerical sclf-consistcnt-field 
technique to find (pop{z) and the energy E{B\^ , p) for each 
plane wave state. 

For a weak or intermediate-strength magnetic field -B|| , 
the effective 2D Hamiltonian takes the form 



2m 
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In Eq.(^, p ~ ~ihX7 — |a(r) is a purely 2D momen- 
tum operator, and p±^ — p- [B\\ x lz\/B\\ is its com- 
ponent perpendicular to _B|| . Two additional terms in 
the free electron dispersion part of H2D are the result 
of the p^-dependent inter-subband mixing. The first of 
them lifts rotational symmetry by causing an anisotropic 
mass enhancement It increases the 2D density of 

states and, for a 2D gas with a fixed sheet density, it 
reduces the Fermi energy calculated from the bottom of 



the 2D conduction band, Ep{B\\) 



7(B|| 



A cu- 



bic term in H2_d is related to the time-reversal symmetry 
breaking by B\\. Note that, depending on the choice of 
a gauge, one may also generate a linear p±^ term, but 
this one can be eliminated by a trivial gauge transfor- 
mation. A perturbation theory analysis of this problem 



is discussed in footnote [Q, and for a moderate field it 
results in the parametric dependences 7 ^ m~^ (Az/As)^ 
and p ^ (Xz/mh) {Xz/Xuf- 

To obtain quantitative estimates for parameters 7 
and /3, we evaluated the electron dispersion at in- 
plane magnetic fields using a full self-consistent nu- 
merical method. The quantum well confining poten- 
tial V{z) was constructed using the nominal growth 
parameters of the sample studied in B which was a 
Al 34Ga eeAs/GaAs heterojunction. V(z) also included 
Hartree and exchange-correlation potentials generated by 
the free carriers in the quantum well. The Hartree po- 
tential was derived from the z-dcpcndcnt 3D density of 
electrons by numerical solution of the Poisson equation. 
The exchange-correlation term was calculated within the 
local-density approximation JT?! . A flat-band boundary 
condition was used, i.e., we assumed that the electric field 
produced by donors in the (Al,Ga)As barrier is screened 
out in the GaAs buffer-layer by the 2D electron gas. 
In each loop of the self-consistent procedure we solved 
numerically the Schrodinger equation with the Hamilto- 
nian in Eq.(^ to get the 3D electron density, neglecting 
M(r, z). Then, a new V{z) was constructed, which entered 
the next loop of the procedure until the self-consistency 
condition was achieved. The numerically obtained de- 
pendences of 7 and (3 on the in-plane magnetic field for an 
electron sheet density of 2 x 10^^ cm~^ are shown in insets 
of Figj^ a) and b). At low fields, 7 ~ and (3 ^ B^^, as 
anticipated in the perturbation theory treatment. The 
proportionality coefficients are plotted in Fig.|l| a) and 
b) versus the electron sheet density. Both the effective 
mass renormalization in the quadratic term of energy dis- 
persions and the time-reversal symmetry breaking cubic 
term are larger at lower 2D electron gas densities, due to 
a weaker confining electric field (i.e., longer A^). 

In H2D in Eq. (^) , disorder is incorporated in the form 
of a scattering potential u[f) « (0|u(r, z)|0). This can 
be characterized by the value of the mean free path, 
I ^ h/pF, or a momentum relaxation time r related 
to the diffusion coefficient D ~ u|r/2. The modifica- 
tion of the electron density of states by i^y only slightly 
affects the value of the electron mean free path. The 
presence of a parallel field also changes the symmetry 
of Born amplitudes of scattering between plane waves 
= e'P^/^^op{z), fpp' = {^Mr,z)\-^p')- Due to the 
momentum-dependent subband mixing, /p^/ acquires an 

fmi^ + ^P^+P'i-W}^ ^here C = 



addition, f^^i 



E„>i glS-p-'gigt-^ -hich is equivalent to the 
presence of a random gauge field in the effective 2D 
Hamiltonian 



a = 2[i3|| X Iz] ^ ■ 

TI>1 



'\u{f, z)|n)(n|z|0) 



The latter can be interpreted as a result of an effec- 
tive 'curving' of a 2D plane by impurities in systems 
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with ^-dependent scattering potential, which in the pres- 
ence of an in-plane magnetic field generates a random 
effective perpendicular field component, bj_ = [rota]^. 
In systems, where scattering is dominated by Coulomb 
centers behind a spacer and is almost independent of 
z, a smaller effect may be taken into account, 6a — 
?7[i?|| x/2]([i?|| X Z^] • V)^u(r). However, (5a has a negligible 
infiuence on the quantum transport characteristics of 2D 
electrons, as compared to the effect of dispersion. 

In the present paper, we study quantum transport 
characteristics of chaotic dots by modelling them as 2D 
disordered billiards with a short-range disorder. It has 
been shown before that the results obtained for a zero- 
dimensional limit of diffusive systems, Tesc ^ L^/D and 
L > I, are universally applicable to the description of WL 
and UCF's in a broad variety of quantum chaotic billiards 
Ml, even in ballistic ones We also used a semi- 

classical diagrammatic language to calculate two-particle 
correlation functions, Cooperons Pc{i^; R, R') and diffu- 
sons P(i{lo; R, R') ||]. These correlation functions emerge 
in the form of ladder diagrams from the perturbation 
theory analysis upon averaging over disorder the Kubo- 
formula conductance. Schematically, the form of a weak 
localisation correction and of the variance and correla- 
tion function of UCF can be represented as gwL(5||) oc 
JdRWiR)Pc{0;R,R) and {{SgiB\\,)Sg{Bp))) cx 
/ dRdR'W{R)W{R')Y.^ c\Pd,c{uj;R,R')\^, where lj = 
-Ef(5||i) — Ef{B\\2), and £'F(i?||) is the Fermi energy 
of the 2D gas calculated from the bottom of the 2D 
conduction band determined in Eq.(^). Dispersionless 
weight factors W{R) both take care of the particle num- 
ber conservation upon diffusion inside a dot p2[ and in- 
corporate coupling parameters to the leads. In the zero- 
dimensional limit, both gwL and {{Sg{B^\i)Sg{Bp))) are 
dominated by the lowest Cooperon (diffuson) relaxation 
mode, + ^C(d) determined by the interplay of the es- 
cape to the reservoirs and the Cooperon (diffuson ) sup- 
pression by time-reversal symmetry breaking (the differ- 
ence in condition of quantum diffusion). The latter effect 
can be analyzed for an infinite system, where the deriva- 
tion is simplified by the use of Fourier representation for 
Cooperons (diffusons). 

Using the Fourier form, the equation for the diffuson 
(Cooperon) correlation function, 11 • Pd(C)(^i'z!) = ''"^^i 
can be obtained from the analysis of a kernel (h — 1), 



dp 

2'KVT 



G\{e,p)GBje-u,±[p-q\), 



where sign '-I-/-' is related to diffuson (Cooperon), respec- 
tively. Disorder-averaged retarded and advanced single- 
particle Green functions, G^'^ correspond to different 
values of B\\. G^'^ were calculated pertubatively with 
respect to all terms containing _B|| , which relies on the 
assumption that within the relevant parametric regime 
the variation of the energy, SE{p-p), induced by is 
small in comparison with the scattering rate, SE ^ h/r. 



The result has the form of the diffusion equation, 
+ S}- DV^ + T-i] Pd = S{R - R'). 



(5) 



It contains ui — uj + S with S — Pf[7(^I1i) ^ 7(^||2)]/2fi- 
and the rate 



-1 ^ 



[7(i?||i) - 7(S||2)]' + —pI[B\\i - S||2]' (6) 



The first term in Eq.(^) comes from the deformation of 
a Fermi circle by B\\ , the second takes into account the 
field effect upon the scattering of plane waves. Eq.(|) also 
contains the difference between the electron kinetic ener- 
gies in two measurements of conductance, w = E-p{B\]^i) — 
Ey{B\\2), each of them shifted, Ey{B\\) = E^-^plj{B\\) 
with respect to the Fermi energy, Ep in the electron gas 
with the same sheet density at B^\ = 0. The latter fact is 
important, since, for lateral dots where electron density 
is fixed, one should substitute uj = uj + S = 0, so that only 
the -dependent anisotropy of the electron wavelength 
along the Fermi line affects the interference pattern of 
current carriers. 

The Cooperon equation derived after the calculation 
of the integral in Tl{uj, q) can be represented in the form 



[-iCj + D{-iV -q)'^ + Tc 1 + T^~i] Pc = S{R - R') 



(7) 



It contains an additional decay rate, ^{B^^i, B|| 



^i_tpI [/3(B||i)-f/3(i?p) 



'B\\ 
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[5||1+S||2]' 



(8) 



which accounts for dephasing between electrons encir- 
cling the same chaotic trajectory in reverse directions, 
the result of lifting the time-reversal symmetry by an 
in-plane magnetic field. This parameter is relevant for 
describing WL correction in Eq. (|l|), and also the o — > u 
crossover in conductance fluctuations. 

To mention, a shift in the Cooperon gauge in Eq.(^, 
q — ^ppm(3{B\\)[B\\ x lz]/B\\, is the result of the fol- 
lowing artifact: cubic term in the effective electron dis- 
persion not only lifts the inversion symmetry of the line 
i?(B||,p) = £^f(S||), but also shifts its geometrical centre 
with respect to the true bottom of the 2D conduction 
band. Since in conductance calculations only electrons 
with E — Ep matter, such a shift would be eliminated 
by choosing a slightly modified initial gauge, which can 
now be corrected by applying a gauge transformation 
Pc — e^^^Pc directly to the Cooperon. One may say 
that, here, the phase-coherent transport is only affected 
by the By-induced {p — > — p)-asymmetric distortion of 
the Fermi circle into an oval, but not by a shift of such 
an oval in the momentum space. 

The effect of time-reversal symmetry breaking by the 
in-plane field described in Eq.(||) has a field dependence 



3 



bB° 



aB?,. In a ballistic billiard, or in a het- 



'B\\ "^\\ .. ^ 

erostructure with z-independent scattering potential, it 
would be dominated by the unusual i?^-dependence (also 
found in ref. [|l3|) that we attribute to the effect of cu- 
bic term generated by in the 2D electron dispersion, 
thus giving rise to a relatively sharp crossover between 
'flat' regions related to orthogonal and unitary symmetry 
regimes. For a large-area (8/im^) quantum dot with elec- 
tron density 2 x 10"cm-2 studied by J.Folk et al ||,|l0|, 
we estimated the crossover field as = 0.6-f0.8T. When 
studying the crossover, one has to take into account that 
the in-plane field also causes fluctuations in conductance, 
without breaking time-reversal symmetry, as described 
by Eqs. For the same parameters of a structure, 

we estimated the field where such a random dependence 
would appear as ~ 0.3T, and the result in Eq. (||) 
suggests that, for a perfectly in-plane field orientation, 
variation of the UCF fingerprint is faster at higher fields. 
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FIG. 1. Calculated dependence of parameters 7 and 
/3 on the sheet density of 2D electrons. Insets show the 
effect of on symmetric, [E{p±) + E{—p±)\/2 and anti- 
symmetric, [E{pi^ — E{—p^_)]/2 parts of the 2D electron 
dispersion in a braoder range of i3| | , where a pertubative 
expansion is not applicable. 
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